Spin Hall effect in a two-dimensional electron gas in the presence of a magnetic field 
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We study the spin Hall effect of a two-dimensional electron gas in the presence of a magnetic field 
and both the Rashba and Dresselhaus spin-orbit interactions. We show that the value of the spin 
Hall conductivity, which is finite only if the Zeeman spin splitting is taken into account, may be 
tuned by varying the ratio of the in-plane and out-of-plane components of the applied magnetic field. 
We identify the origin of this behavior with the different role played by the interplay of spin-orbit 
and Zeeman couplings for in-plane and out-of-plane magnetic field components. 



The spin Hall effect is the transverse pure spin current 
response to an applied electric field, with the spin current 
polarization perpendicular to both the spin current and 
the electric field. This effect in semiconductors is due to 
the spin-orbit coupling: for a review see Ref. [l[ . In a dis- 
ordered two-dimensional electron gas (2DEG), it is now 
accepted that, at least in uniform and stationary con- 
ditions, the spin Hall conductivity (SHC) vanishes for 
the Rashba spin-orbit interaction or linear Dresselhaus 
spin-orbit coupling 0, 0, Gl> HI- The reason for this has 
been clearly associated to the linear-in-momentum de- 
pendence of the spin-orbit interaction in a 2DEG 0,0,1. 
As pointed out by Rashba 0, disorder is only important 
in order to obtain a steady state and the vanishing of the 
SHC persists in the presence of a perpendicular magnetic 
field if the Zeeman splitting is neglected. On the other 
hand, in a disordered 2DEG, the presence of the Zee- 
man coupling due to an in-plane applied magnetic field, 
may lead to a non- vanishing SHC [J 1(J HI- Therefore, 
it is relevant to further investigate how the interplay of 
Zeeman coupling and various spin orbit interactions may 
lead to a finite SHC in a 2DEG in a magnetic field even 
in the absence of disorder. In this paper we show that the 
role played by the Zeeman coupling depends on whether 
the applied magnetic field has an in-plane component in 
addition to the one perpendicular to the plane of the 
2DEG. 

In the case of a perpendicular magnetic field, Rashba 
provided an elegant proof of the vanishing of the SHC 
by deriving a sum rule based on the knowledge of the 
exact eigenvalues and eigenstates of the Hamiltonian [3] • 
When both the Rashba and Dresselhaus interactions are 
present an exact solution for the eigenstates and eigen- 
values is only known in terms of recurrence relations for 
the expansion coefficients in the Landau states basis [3]. 
For the special case when the two spin-orbit interactions 
have equal strenghts one can again explicitly diagonalize 
the Hamiltonian 14 , 1^, 1(| 17 1 . In this case the SHC has 
been investigated in the presence of Zeeman coupling (l8|. 
However, when the Zeeman coupling is due to a perpen- 
dicular magnetic field, the SHC is associated to a charge 
Hall conductivity (CHC). As the Zeeman coupling gives 



rise to a perpendicular spin polarization, it is not surpris- 
ing that the CHC also carries a spin current. In the case 
of a tilted magnetic field, we show that, in a clean sample, 
a new contribution to the SHC appears and, most inter- 
estingly, it is anisotropic with respect to the direction 
of the in-plane component. This phenomenon is similar 
to what has been predicted for the case of a disordered 
2DEG. 

The Hamiltonian for the 2DEG reads 

H = H Q + H Z + H so , (1) 

where the spin-independent part is given by 
1 



H = —IT, n = p + -A, e > . 

2m c 



(2) 



Here the electrons are supposed to be confined in the 
x — y plane. It is useful to choose, for an arbitrarily 
oriented magnetic field B = (B x ,B y , B z ), the gauge A = 
(—B z y + B y z, —B x z, 0). The Zeeman coupling is 



H z = ^g^B'B ■ a-, 



(3) 



where er = (a x , <7 y , a z ) is the Pauli matrix vector. Fi- 
nally, the spin-orbit Hamiltonian is given by 

Hso = -T {H-yCTx - II^CTy) + - (HyCTy - T1 x <j x ) , (4) 

with the two terms describing the Rashba and the linear 
Dresselhaus spin-orbit interactions, respectively. 

Due to our gauge choice, the system has translational 
invariance along the x-direction, p x = hk is a good quan- 
tum number and the orbital part of the eigenfunctions 
can be chosen as 



x,y) = e lkx ip(y-y ), 



(5) 



where ip is a generic harmonic oscillator wavefunction 
and 7/o = ' 2 fc indicates the orbit center coordinate with 
lb = y/hc/eB the magnetic length. By introducing the 
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creation and annihilation operators for the harmonic os- 
cillator 



Tl x + Hl y = -\f2M b a) 



iH y = —^/2hlt,a : 



(6) 



and by measuring all the energies in units of the ciclotron 
energy Huj c — heB z /mc, we get the Hamiltonian |T]) in 
the following dimensionless form 



H 



(at a + 1/2) + f ( 



B x 
K 



El 

B, 



(J) 



R [a 1 a 



V2r] D (a)a + 



where g R = amlb/ti 2 and t]d = Pinl^/h 2 . In order to 
analyze the role played by the various components of the 
magnetic field, it is useful to introduce g x — gB x /B z , 
g y = gB y /B z , and g z = g to be treated as independent 
parameters. 

Following Rashba we compute the zero frequency SHC, 
whose Kubo formula reads 



-— -j hm 

b A, A' 



3m(<AK|A'>(A'|J M |A)) 



^A 



(8) 



where |A) are the eigenstates with eigenvalues E\ — huj\ 
and the star above the sum indicates that we are only 
considering states |A') below the Fermi energy (here set 
to zero for simplicity) and |A) above the Fermi energy. 

The spin current operator, for the z-component of the 
spin, is defined as 



Jzi = -^{vzVi + v l a z ). 



(9) 



In terms of the harmonic oscillator operators of Eq.([6]), 
the velocity operators can be derived according to v x — 
[x, H]/ih and v y — [y, H]/ih. They read 



h 



iy/2h 



(a f +a + V2g R a y + V2rj D a x ) (10) 
(i(at - a) + V2rj R a x + V2g D a v ) , 



so that the spin current becomes 



Jzx — 



iV2l b 



+ a) a z . 



(11) 



Eq.© together with Eq. (flU|) is a generalization of 
Eq.(5) of 0| . In Ref.Q the author studies a case with 
only the Rashba spin-orbit interaction and a magnetic 
field orthogonal to the 2DEG in the absence of Zee- 
man splitting. In that case, the matrices (AI^IA') and 
(A' | J zx | A) are antisymmetric and symmetric, respectively 
(with respect to the exchange of A and A'), giving rise to 
a purely imaginary matrix element (A|u 3/ |A')(A'| J zx \ty- 
Hence the two terms in the commutator appearing in the 



Kubo formula are equal and opposite giving just a fac- 
tor of two. By contrast, we study the more general case 
in which also the Dresselhaus and Zeeman couplings are 
present. In this case no general argument about the sym- 
metry of the matrix elements holds and, therefore, the 
imaginary part in Eq.© has to be taken explicitly. 

In order to numerically diagonalize the Hamiltonian we 
expand the eigenstates |A) = J2 n <r c no-l ncr ) m the Landau 
levels (LL) basis 



(r) = -^tMO) <g> \cr), 



(12) 



corresponding to the eigenvalues E n = huj c (n + 1/2) in 
the absence of spin-orbit interactions and Zeeman cou- 
plings. By working out the matrix elements, Eq.(|8|) for 
the SHC becomes 



1 eui, 



■o 



1 



2 nh fv ( WA ' 



Re(4u, * B xx ,), (13) 



E 



+ ^+Tctc^ +la , (14) 



n' 

- ^^ 7 +lc A ;; li c A u - Vrttf+n&J. (15) 



^zxy (0) 

where 
B\\> 



The numerical evaluation of Eq. (fT3| shows that, in the 
absence of the Zeeman coupling, the SHC vanishes in all 
cases: whether the Rashba, the Dresselhaus or both spin- 
orbit couplings are present. This is a generalization of 
what obtained in Ref.0]. In that case it has been shown 
that the SHC is zero in the presence of Rashba spin-orbit 
coupling, with a sum rule argument. That result is also 
valid in the case of a Dresselhaus spin-orbit coupling, as 
this kind of interaction can be obtained, starting from 
the Rashba one, by use of the following rotation in the 
spin pace 



(Tj => <Xy((7~ => — i<J + ) 
(7, => (T z . 



(16) 



However, in spite of the fact that, when both interactions 
are present, we were not able to extend the sum-rule ar- 
guments, our numerical calculation clearly shows that 
also in the presence of both kinds of spin orbit interac- 
tions the SHC vanishes. 

In order to study the effect of the Zeeman coupling on 
the SHC, we consider first the case with the magnetic 
field orthogonal to the 2DEG with no in-plane compo- 
nent. In Fig[T] we plot the low lying part of the energy 
spectrum (upper panel) and the SHC (lower panel) as a 
function of g z . The energy spectrum shown corresponds 
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^R= U % 
•T1 R =0.1T1 D =0 

,t! r =0.2ti d =0 
•T1 R =0 r, D =0.1 
,r, R =0 T! D =0.2 



FIG. 1: (color on-line) (Upper panel) Low lying energy spec- 
trum as function of the Zeeman coupling g z . The arrows in- 
dicate the spin polarization of the LL's. Red arrows indicate 
the unbalanced spins below the Fermi levels. (Lower panel) 
SHC in units of e/Anh as function of g z for different values of 
r}n and r\o (indicated in the plot). lo c = 2meV 



to the case r/R = = 0. For simplicity, we have chosen, 
at g z = 0, a Fermi level lying between the second and the 
third LL, (the dashed line in the figure). Here, and in the 
following, the cyclotron frequency is uj c = 2meV corre- 
sponding to a magnetic field of ~ IT for GaAs based 
materials. The corresponding SHC for g z = is given by 
the black curve in the lower panel. It is quantized be- 
cause the CHC is, due to the choice of having an integer 
number of LLs fully occupied. While the CHC is constant 
with g z , the SHC shows sharp steps at those level cross- 
ings at which the total spin polarization of the Fermi sea 
changes. To emphasize how this happens we have con- 
nected in FigfTJby vertical lines the energy level crossings 
with the steps of the SHC. We can clearly see, by inspec- 
tion, that only an odd number of crossings gives rise to 
a new net polarization of the current carrying states, re- 
sulting in the steps of the SHC of Fig[TJ Of course, here 
spin and charge conductance quantizations go together. 
The other curves in the lower panel of FigQ] describe the 
behavior of the SHC at finite values of tjr [19( or tjd and 
have been shifted uniformly for clarity. A finite Rashba or 
Dresselhaus slightly modifies the shape of the SHC steps, 
without substantially changing the scenario obtained in 
the absence of spin-orbit interaction. In particular, there 
is a rounding effect and a small shift of the g z values at 
which the steps occur, which is positive or negative in the 
case of Dresselhaus and Rashba coupling, respectively. 

Next we turn to the case in which the magnetic field 
is tilted with respect to the z-axis and has both in-plane 
and out-of-plane components. Two different subcases are 
considered: a magnetic field lying in the xz plane or in 



the yz plane. Moreover, in order to emphasize the effect 
of the in-plane Zeeman coupling, we will first neglect the 
effect of g z . Therefore the relevant Zeeman parameters 
are g x and g y . These cases are analysed in detail in Figs. 

1E1S1 
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b-b T1 R =0.08 r, D =0.2 
-^r| R =0.16r, D =0.2 




FIG. 2: (color on-line)Spin Hall conductivity in units of e/Anh 
as function of gy or gx for different values of gn and t]d 
(indicated in the plot). 



In Fig[2]we plot the SHC at a fixed r\R ^ for various 
r\D (lower red curves) and at a fixed ^ for various 
(top black curves). The former are plotted as a function 
of g y and the latter as a function of g x . This is because at 
r]B = 0, i]d 7^ the SHC is identically zero if g x = for 
any g y ^ 0, while at tjr ^ 0, r\v — 0, the SHC vanishes 
for any g x ^ if g v = 0. The vanishing of the SHC at 
9y — 9x = has been found in ref.[7|. At increasing g y { x ) 
the peak in the SHC is larger the larger is the difference 
between 77^ and t\b- 




FIG. 3: (color on-line)SHC as a function of H)r and go for 

g x = 0, g y = 2 
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FIG. 4: (color on-line)SHC as a function of rjn and t/d for 

g x = 2, g y = 




The curves in Fig|2] show a remarkable reflection sym- 
metry with respect to the x-axis. Indeed, the Rashba 
Hamiltonian, in the presence of a g y Zeeman splitting, 
maps into the Dresselhaus Hamiltonian in the presence 
of a g x Zeeman splitting. However this property is even 
more general, provided that g z vanishes, as the unitary 
transformation of Eq. (|16[) maps the Hamiltonian into 
itself by interchanging the pair (r}R,g y ) with the pair 
(VD,gx)- While the velocity operator v does not change 
under this transformation, the current J zy changes its 
sign. As a result the SHC is antisymmetric with respect 
to the interchange of rjn,g y with rjD,g x and therefore 

^zxy (VR = a > Vd = b, g x = c, g y = d) = (17) 
= - ^zxy (r)R = b, i]d = a, g x = d, g y = c). 

This result has already been found in the absence of mag- 
netic field In Figs. [3l[4l the color maps refer to the 
SHC for a wide window of t]r, rju (x and y axis). In Fig[3] 
g x — and g y — 2, while FigQ] shows the complementary 
case g x = 2 and g y = 0. The color pattern in the figures 
[3] and 5] is exactly the same but with the interchange of 
the sign of the SHC because of the symmetry property 
(fP7|) . It can be noticed that the SHC can be sizeably 
different from zero in a wide range of parameters. We 
stress that the SHC shown here is not associated with 
a spin polarized charge Hall current as, up to now, we 
have been considering just the case g z = 0. As g x and g y 
include the components of the magnetic field in their def- 
inition, exchanging g x and g y is equivalent to a change of 
direction of the in-plane magnetic field. This implies that 
when the magnetic field is in the x direction a nonzero 
SHC only occurs for finite 77/5, while when the field is in 
the y direction the SHC occurs for finite r/R only. 

However, in a real experiment of a 2DEG with both 
perpendicular and parallel magnetic fields, one has also 
the Zeeman spin splitting in the z-direction parametrized 
by a non- vanishing g z . We conclude, then, our analysis, 
by discussing this case. We choose the magnetic field 



FIG. 5: (color on-line) Non additivity of the SHC with g y and 



oriented in the y — z plane, as the x — z orientation can 
be inferred by symmetry, according to Eq. (|17p . The ref- 
erence curve, for g z = is plotted in the lower part od 
FigGH The SHC is non vanishing only in a restricted 
region of g y values because, in the absence of g z , the 
electron spins tend to polarize along the y-axis, when g y 
grows, so inhibiting the Rashba spin precession. As a 
result the SHC tends to zero since no z-axis polariza- 
tion is available for the spin current. This does not hap- 
pen when g z is finite as it appears in the upper curve of 
FigEl there, the competition of the two orthogonal Zee- 
man splittings {g y and g z ) does not allow for an in-plane 
polarization effect, therefore letting the Rashba (or Dres- 
selhaus) coupling playing its role. The measured SHC, 
in this case, will have two contributions: the first one is 
due the Zeeman polarization of the LL discussed so far 
in Fig. [TJ the second (the most interesting one) is the 
SHE due to the spin-orbit coupling of Figs l2l3l4l As it 
can be seen from the Fig. O the two effects do not su- 
perimpose trivially, i.e., the SHC in the presence of both 
a Zeeman interaction in the z and y direction is not the 
linear superposition of the SHC with g y or g z only. To 
stress this point we have plotted the SHC as a function 
of g y in Fig. with g z = (black curve) and g z = 2 
(red curve) for i]r = 0.2 (and i]o = for simplicity). 
The behavior of the two curves as function of g y is sub- 
stantially different. The striking result is that a finite g z 
makes the SHC finite even at increasing values of g y , in 
a region where the SHC would vanish if g z = 0. This is 
clearly seen in the inset of Fig. 03 where the red curve 
has been shifted downward to eliminate the background 
SHC associated with the cohexisting CHC with a finite 
g z . Of course, the derealization of the orientation of the 
spin, due to presence of both non commuting spin terms 
in the Hamiltonian, makes the tilting of the spin in the 
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z-direction easier even at larger g y . 

In a real sample, one expects that both g y (x) and g z 
are different from zero and in a nonzero SHC can be 
expected when a tilted magnetic held is applied. 

In conclusion, we have shown that a hnite spin Hall 
effect is possible when both out-of-plane and in-plane 
magnetic field components are applied to a 2DEG in the 
presence of Rashba and/or Dresselhaus spin-orbit cou- 
plings. The interplay of the Zeeman couplings along the 
perpendicular and parallel direction, with respect to the 
plane of the 2DEG, allows for a tuning of the SHC and 
should provide the way for an experimental detection. 

Depending on whether the Rashba or Dresselhaus spin 
orbit coupling dominates, the SHC changes with chang- 
ing of the orientation of the in-plane component of the 
magnetic field. 
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